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^ ■ Abstract 

1^^ I In this paper we consider a representative a priori unstable Hamiltonian system 

^SJ ' with 2 + 1/2 degrees of freedom, to which we apply the geometric mechanism for 

[~^ . diffusion introduced in the paper Delshams et al., Mem. Amer. Math. Soc. 2006, 

^^ I and generalized in Delshams and Huguet, Nonlinearity 2009, and provide explicit, 

concrete and easily verifiable conditions for the existence of diffusing orbits. 

The simplification of the hypotheses allows us to perform explicitly the computa- 
tions along the proof, which contribute to present in an easily understandable way the 
^ ' geometric mechanism of diffusion. In particular, we fully describe the construction 

j^ ' of the scattering map and the combination of two types of dynamics on a normally 



hyperbolic invariant manifold. 



1 Introduction 



The goal of this paper is to apply the geometric mechanism for diffusion introduced in 
|DLS06j and generalized in |DH09j , to a representative a priori unstable Hamiltonian 
system with 2 + 1/2 degrees of freedom, establishing explicit conditions for diffusion. 



* Supported in part by MICINN-FEDER MTM2009-06973 and CUR-DIUE 2009SGR859 grants. 



The phenomenon of global instability in nearly integrable Hamiltonian systems, 
called Arnold diffusion, deals essentially with the question of what is the effect on the 
dynamics when an autonomous integrable mechanical system is subject to a small 
periodic perturbation. That is, whether these effects accumulate over time leading 
to a large term effect or whether they average out. 

For an integrable Hamiltonian system written in action-angle variables all the 
trajectories lie on invariant tori, with an associated dynamics consisting of a rigid 
rotation with constant frequency. For a general perturbation of size e of a non- 
degenerate integrable Hamiltonian, the KAM Theorem (see |Lla01] for a survey) 
ensures the stability for most of the trajectories of the system. More precisely, those 
invariant tori in the unperturbed system e = having Diophantine frequencies are 
preserved (they are tori with non-resonant frequencies), giving rise to a Cantorian 
foliation of invariant tori for the perturbed system e > 0. Thus, the set of surviving 
tori has a large measure but also many gaps among them, which can be of size up 
to order ^/£. However, nothing is said by the KAM theorem about the stability 
of the trajectories that do not lie on the non-resonant invariant tori. Besides, for 
systems with more than two degrees of freedom the invariant tori are not anymore 
an obstruction for the existence of trajectories that may possibly drift arbitrarily far, 
called diffusing orbits. 

The first description of this instability phenomenon was given by Arnold in 
|Arn64j by means of his celebrated example and has been thoroughly studied since 
then, using a wide range of techniques: geometric, variational and topological (see 
[DGLSOS] for a long list of references) . 

The geometric mechanism for diffusion in |DLS061 IDH09J is based on the classi- 
cal Arnold mechanism for diffusion, which consists of constructing a transition chain, 
that is, a finite sequence of whiskered transition tori (lower dimensional invariant tori 
having non-resonant frequencies with associated stable and unstable manifolds) with 
transverse heteroclinic intersections, that is, the unstable manifold of each transition 
torus intersects transversally the stable manifold of the next one. Arnold |Arn64j 
considered an integrable system with a hyperbolic component (a rotor and a pendu- 
lum) plus a particular periodic in time perturbation, that allowed him to construct 
a transition chain of primary KAM tori, that is, lower dimensional tori that are 
just a continuation of the lower dimensional invariant tori that existed in the unper- 
turbed case. Using a topological argument, he proved that diffusing orbits exist in a 
neighborhood of the tori in a transition chain. 

Nevertheless, when one considers a generic perturbation, one faces the problem 
that the gaps of size ^/e in the foliation of primary KAM tori are bigger than the size e 
of the heteroclinic intersection of their whiskers. Therefore, this prevents the unstable 
whisker of a primary KAM torus of this foliation intersecting the stable whisker of the 
next surviving primary KAM torus and one can not construct a transition chain using 
only primary KAM tori. This is known in the literature as the large gap problem and 
has been solved very recently by different methods |DLSnni IDLSOBl IDHn9[ l(:;Yn4[ 
[CY091 nviOil IPTOTI iGLOebl IC;L06a| . 

The strategy in [DLSOG) IDH09J to overcome the large gap problem was to incor- 
porate in the transition chain other invariant objects, which are not present in the 



unperturbed system and are created by the resonances therein, in order to fill the 
gaps between two primary KAM tori. Among them, there are the so-called secondary 
KAM tori, which are lower invariant KAM tori topologically different from the pri- 
mary ones. The scattering map jDLSOSj is the essential tool to study the heteroclinic 
connections between invariant objects like primary or secondary KAM tori. 

In |DLS06j it was proved the existence of Arnold diffusion in a priori unstable 
Hamiltonian systems of 2 -|- 1/2 degrees of freedom, under concrete geometric hy- 
potheses. However, one of them was the assumption of a non-generic condition, 
namely, that the Hamiltonian was a trigonometric polynomial in the angular vari- 
ables. This latter assumption was removed in |DH09] and therefore, for these kind 
of Hamiltonians, the conditions required for the geometric mechanism of diffusion 
were proven to be generic in the C^ topology. Moreover, the removal of this hypoth- 
esis allowed us to present the conditions for diffusion more explicitly in terms of the 
original perturbation. 

A strong feature of this mechanism, in contrast to other existing ones, is that the 
conditions for diffusion are computable and therefore verifiable in specific examples. 
Moreover, the way the mechanism is designed allows us to give an explicit description 
of the diffusing orbits. From our point of view, this fact makes this mechanism really 
attracting for applications, where the computation of the diffusing orbit in concrete 
systems is the cornerstone of the problem (see JDMR08] ). 

Although the conditions for diffusion for any concrete system are totally explicit 
and computable, the computations needed involve the application of several steps 
of the averaging method, the expansion in e of a NHIM, and the verification of the 
existence of non-degenerate critical points of the Melnikov potential along the straight 
lines. The computation and verification of these conditions may hide the elementary 
features of the geometrical method. Because of this, in this paper we have chosen 
a representative kind of a priori unstable Hamiltonian systems, where the required 
hypotheses for the application of the geometric method are trivially fulfilled. On the 
contrary, the geometry of the existence of non-degenerate homoclinic orbits to the 
NHIM as well as the behavior of the scattering map on the NHIM and its interaction 
with the inner dynamics in the NHIM can be fully described. 

The main result of this paper is Theorem 12.11 which states the existence of diffu- 
sion under very concise and easily verifiable hypotheses, for a representative class of 
priori unstable Hamiltonian systems with 2 + 1/2 degrees of freedom. 

Our aim is to explain the strategy of the mechanism in a clear and understandable 
way, and for this reason we have chosen an illustrative type of Hamiltonian systems 
which appear commonly in the literature, for which the computations along the proof 
can be performed explicitly. In Section [2] we discuss the reasons for the particular 
choices we have made. Overall, we wanted to avoid the problem that the technical 
details hide the main ideas behind the mechanism. 

Another important point of this paper is that we can give explicit expressions for 
the equations defining the invariant tori as well as for the Melnikov function and the 
reduced Poincare function, which are essential for the computation of the scattering 
map. For those readers interested in further their understanding of the scattering 
map, will find here a good illustration of its construction for a particular example. 



based on geometrical considerations. 

Moreover, in our description we try to present an approach that emphasizes those 
points that were crucial to prove diffusion. Hence, we describe the two different 
dynamics on a normally hyperbolic invariant manifold that need to be combined to 
create diffusion and we compute them explicitly. 

Although this paper strongly relies on the results obtained in the previous papers 
[DLS061 IDH09J . we have made an effort to make it self-contained for the reader 
mainly interested in the heuristic description of the mechanism and how it applies 
to concrete examples. However, we accompany the exposition with precise references 
to the results in |DLS061 IDH09J for those readers concerned with the rigorous proofs 
for more general systems. 

The paper is organized as follows: in Section [2] we introduce a representative 
kind of a priori unstable Hamiltonians with 2 + 1/2 degrees of freedom considered in 
this paper, for which we can state our main result. Theorem 12.11 which establishes 
the existence of a diffusing orbit for the model considered. In Section [3] we perform 
the explicit verification of the mechanism for the Hamiltonian of Theorem 12.11 The 
verification is structured in four parts, and includes a detailed description of the 
scattering map. 



2 Set up and main result 



In |DH09j there were given explicit conditions for the existence of a diffusing orbit 
for generic a priori unstable Hamiltonian systems. That paper was mainly focused 
on proving the genericity of the result, so although conditions were explicit, some 
computational effort was required to check them. As we already mentioned in the 
introduction, in this paper we plan to give sufficient conditions, easily verifiable, 
which guarantee the existence of diffusion for a general case of a priori unstable 
Hamiltonian systems. 

In this section, we first introduce a representative type of a priori unstable Hamil- 
tonian systems of 2 + 1/2 degrees of freedom, which is usually found with several vari- 
ations in the literature |(]hi79[ S7.5] |(X;94[ §9, S12] |B(:;vni[ IBBn2[ iBBBOHl Il'ren4j . 
when explicit computations are performed. We first discuss the features and partic- 
ularities of this type of systems and we finally state Theorem 12. H which establishes 
the existence of diffusing orbits for these systems under explicit and easily verifiable 
conditions. 

We consider an a priori unstable Hamiltonian system as introduced by Chier- 
chia and Gallavotti |CG941 sections 7.5 and 12], consisting of a 27r-periodic in time 
perturbation of a pendulum and a rotor described by the following non-autonomous 
Hamiltonian 

He{p,qJ,y^,t) = Ho{p,q,I) + eh{p,q,I,ip,t;e) 

1 (1) 

= P±{p,q) + -I^ + eh{p,q,I,ip,t;e). 

We notice that a motivation for the model above comes from a normal form around a 
resonance of a nearly integrable Hamiltonian, and we refer the reader to |DG0HlKL08j 
for more details. 



The second term i^l"^ in the integrable Hamiltonian Hq{p, q, I) of Hamiltonian ([T]) 
describes a rotor and the first one 

P±iP,q) = ±(^\p' + Viq)^ (2) 

a pendulum, where the potential V{q) is a 27r-periodic function, whose non-degenerate 
maxima give rise to saddle points of the pendulum ^ and therefore, to hyperbolic 
invariant tori of the Hamiltonian ([T|) when e = 0. Typically it is assumed that the 
maximum of V is attained at the origin q = 0, where V is assumed to vanish, as it is 
case for the standard pendulum, where 

V{q)=cosq-l. (3) 

This is the simple and standard choice of potential V{q) that we are going to consider 
in this paper, so that 

P± {p, q) = ± (pV2 + cos g - 1) . 

The origin (p = 0, g = 0) is a saddle point of the standard pendulum, and its 
separatrix P^ (0) for positive p is given by 

qo{t) = 4arctane'^*, po{t) = 2/cosht. (4) 

Notice that other choices of V give rise to different separatrices that are not usually 
so simple. 

The negative sign in the pendulum ([2]) is only to emphasize the fact that the 
geometric mechanism we are using does not require the Hamiltonian Hq to be positive 
definite, as it is the case in the variational approach, see |CY09j . 

The term eh in ([TJ is the perturbation term and depends periodically on time 
and on the angular variable ip, so that h can be expressed via its Fourier series in the 
variables {(p,t) 

h{p,q,I,y,,t;e)= J^ h,^i{p,q, I-e)e"^''^+''\ (5) 

It is common in the literature to consider a perturbation ([5]) depending only on the 
angular variables {q,ip,t), and formed by the product of a function of the pendulum 
variable q times a function of the angular variables {(p, t) 

Hp, q, I, 'P, t) = f{q)9{^, t)- (6) 

This is the kind of perturbation that we are going to consider along this paper, 
particularly because with this choice of h the Melnikov potential (jlSp , which will be 
an essential tool for the computations along the paper, has the same harmonics as 
the function g, and they can be computed explicitly. So we will focus on a concrete 
type of Hamiltonians of the form 

H,{p,q,I,^,t)=±i^ + cosq-l\+ — + ef{q)g{ip, t), (7) 



defined for any real value of {p, q, I, 93, t, e) and 27r-periodic in the angular variables 
{q,ip,t). 

The function / could be any 27r-periodic function. However, to easily compute the 
harmonics of the Melnikov potential ()15|) . we are going to assume along this paper 
that / has the simple form: 

f{q)= cos q. (8) 

About the choice of / we would like to remark two important features. On the 
one hand, thanks to the fact that /'(O) = the problem does not require the use 
of the theory of normally hyperbolic invariant manifolds to ensure the persistence of 
these type of objects. This simplifies the exposition and the computations, but since 
we do not assume /(O) = the problem maintains all the richness and complexity of 
the large gaps problem. So, although the choice dSD for / may seem very restrictive, 
we would like to insist on the fact that the complexity of the original problem is 
preserved. At the beginning of Section [XT] we discuss in detail the role of the function 
/ in the problem. 

A general function g is of the form 

9{^,t) = X] ak,icos{k^-lt-ak,i) + dk,icos{kip + lt- ak,i), 

(fc,0eN2 

which in general has an infinite number of harmonics in the angles (if, t) and where 

Since for simplicity we will study diffusion only for / positive along this paper, 
we will consider ak,i = 0, that is, 

g{ip,t) = ^ ak,icos{lnp-lt-ak,i). (9) 

(fe,0eN2 

In a natural way, and also for simplicity, we have chosen g to be an analytic 
function and therefore we will assume an exponential decay for its Fourier coefficients. 
That is, \ak^i\ < e~ '^^''''', where 6 is the size of the domain of analyticity. In this 
paper we simply are going to assume that they have some general lower bound with 
exponential decay, that is 

where 1 < /? < 2. Or, equivalently, we are going to assume 

q/V'^' < |ofc,/| < o/r', (10) 

where 1 < /3 < 2 and < a < a. Moreover, < p,r < 1 are real numbers that will 
be chosen small enough. 

Although the lower bound for the coefficients ai^i in the above equation seems 
very restrictive, there are several reasons for this particular choice. For the more 
expert reader, let us say that condition (jlOp implies that big gaps of maximal size 
appear for all the resonances in first order with respect to the parameter e, without 
performing any step of averaging. This feature is explained thoroughly in Section 
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after equation (j42p . A second reason is that requirements (jlOp are simple to state 
and verify. A generic, and, of course, more technical, set of conditions for generic 
perturbations are given explicitly in |DH09j . When the lower bound condition for a^i 
in (fTO|) are not satisfied, several steps of averaging are needed. 
We can now state our main result: 

Theorem 2.1 Consider a Hamiltonian of the form ^, where f{q) is given by ^ 
and g{ip, t) is any analytic function of the form ([9]) with non-vanishing Fourier coef- 
ficients satisfying (jlOp . Assume that 



A 



ai,o 



ao,i 



< 1/1.6 or A > 1.6, (11) 



as well as < p < p* and < r < r* , where p*{X,a,a, 13) and r*(A,a,a,/3) are small 
enough. 

Then, for any I^ > 0, there exists e* = £*(/+) > such that for any —1/2 < 
/_</+< /ijl and < e < e* , there exists a trajectory {p{t),q{t),I{t),(p{t)) of the 
Hamiltonian ^ such that for some T > 

/(O) < /_; I{T) > /+. 

We want to remark now that not every perturbation ([6|) gives rise to diffusion 
in the action I. In particular, if the function g{ip, t) in ([9]) does not depend on ip, 
the action / is a first integral, so it does not change at all. Moreover, if g{(p, t) does 
not depend on t, Hamiltonian ([7|) is autonomous and therefore H^ is a first integral, 
so that only deviations of size ^/£ are possible for the action /. The same happens 
when the function g{ip, t) does not depend fully on the two angular variables, but 
only through a integer linear combination of them — that is, g{(p,t) = G{ip), where 
ip = kQ(p — lot is an integer combination of the angular variables (/?, t) — , as can be 
easily checked by introducing ^ as a new angular variable. In these three cases, an 
infinite number of Fourier coefficients a^^i of the function g{ip, t) in ([UJ vanish. 

Remark 2.2 The condition (|lip for ai.o and ao.i jointly with the assumption p*,r* 
being small enough, will ensure the existence of a global connected homoclinic mani- 
fold in Section \3.1[ Indeed, for p* , r* small enough, thanks to the upper bound (|lUp , 
the perturbation g in Q can be approximated by its first order trigonometric polyno- 
mial 

gr[<i] (y,^ t) = ao,o + ai,o cos((/? - cJi^o) + ^0,1 cos(s - o-o,i), 



in such a way that all the computations required for the Melnikov potential (J15p will 
depend explicitly on g^-^'. 



3 Proof of Theorem 12.1 



In the following sections, we will consider any Hamiltonian satisfying the hypotheses 
of Theorem 12.11 and we will show how the geometric mechanism in |DLS061 IDH09J 
can be applied to construct diffusing orbits. 



We already mentioned in the Introduction of this paper that the geometric mech- 
anism in |DLS061 IDH09J is based on the classical Arnold mechanism for diffusion. 
It consists of constructing a transition chain, that is, a finite sequence of whiskered 
transition tori such that the unstable manifold of each torus intersects transver- 
sally the stable manifold of the next one. As a main novelty, in |DLS06l IDHOQj the 
transition chain incorporates primary as well as secondary KAM tori created by the 
resonances (as already mentioned in the introduction), in order to overcome the large 
gap problem. 

In this paper, we will try to present a description of the mechanism that empha- 
sizes more the geometrical aspects. We think that this description may contribute to 
a better understanding and applicability of the mechanism. 

In order to prove the existence of a diffusing orbit we will identify first a NHIM 
(normally hyperbolic invariant manifold) with associated stable and unstable mani- 
folds). It will organize the different invariant objects involved in the transition chain 
(the skeleton for the diffusing orbit). 

The diffusing orbit we are looking for starts on a point close to the NHIM and 
in finite time reaches another point close to the NHIM but arbitrarily far from the 
original one. Of course, if the starting point lies just on the NHIM (3D), the invariant 
tori (2D) inside the NHIM act as barriers for diffusion and the orbit is confined in 
a bounded domain. Fortunately, there exists an external dynamics to the NHIM, 
provided by its associated stable and unstable manifolds, which will be essential to 
overcome the obstacles of the invariant tori and escape from them, as long as the 
starting point does not lie on the NHIM but very close to it. Hence, it is crucial 
for the mechanism that the external dynamics does not preserve the invariant tori 
existing in the NHIM. Otherwise, the orbit will be confined in a finite domain by 
both the inner and the external dynamics with no possibility to escape. 

Therefore, given a Hamiltonian of the form ([7]), the proof of the existence of 
diffusing orbits can be sketched with the following steps: Detect the NHIM (3D) and 
the associated stable and unstable manifolds (4D), determine the inner and the outer 
dynamics of the NHIM as well as the invariant objects for each one, and finally show 
that the outer dynamics does not preserve the invariant objects for the inner one. 

One of the novelties of this paper is the explicit description of the outer dynamics 
provided by the scattering map [DLSOSj . It is given by the time-e flow of a Hamil- 
tonian that in first order is given by an autonomous Hamiltonian of one degree of 
freedom, therefore integrable. Moreover, using a geometric description, we are able to 
obtain an explicit expression for this autonomous Hamiltonian, which is the reduced 
Poincare function ([35|1 with the opposite sign. 

On the other hand, using averaging theory, one can show that the Hamiltonian 
defining the inner dynamics can be transformed into a normal form consisting of an 
integrable Hamiltonian plus a small perturbation. 

Thus, we have two dynamics defined on the NHIM that can be approximated in 
suitable coordinates by one-degree of freedom autonomous Hamiltonians. Therefore, 
the invariant objects are given approximately by the level sets of these integrable 
Hamiltonians, for which we provide explicit expressions. 

Finally, we impose that the scattering map moves the invariant tori for the inner 



dynamics, in such a way that the image under the scattering map of each of these 
invariant tori intersects transversally another torus. We will show that this is a 
generic property. 

3.1 Part 1. Existence of a NHIM and associated stable 
and unstable manifolds 

The first part of the proof consists of looking for a NHIM with associated stable and 
unstable manifolds that intersect transversally. In order to prove the existence of 
these invariant objects, we will look for them in the unperturbed case, which is much 
simpler, and then study the persistence of these objects under the perturbation. 

For e = 0, Hamiltonian Hq in d?]) consists of two uncoupled systems: a pendulum 
plus a rotor. Therefore, it is clear that the 3-dimensional manifold given by 

A = {(0, 0, /, if, s) : (I, ip,s) £Rx T^} (12) 

is an invariant manifold with associated stable and unstable manifolds (inherited 
from the separatrices of the pendulum). These manifolds coincide along a separatrix 
given by 

W'A = W^A = {{po{T),qo{T),I,ip,s) :T€R,Ie [-1/2,/;], (v9,s) G T^}, (13) 

where {po{t)jQo{t)) is the chosen orbit ^ of the pendulum, which is homoclinic to 
the saddle point p = 0,q = 0. 

The integrable Hamiltonian Hq has a one-parameter family of two-dimensional 
whiskered tori given by 

7? = {(0,0,/,v9,s):(^,s)gT2}, (14) 

with associated frequency (1, 1). 

When we consider the perturbation h, that is e > 0, using the standard theory of 
NHIM, see |HPS771IFen79j . we know that for e > small enough, maybe restricting 
A to a compact subset, the manifold A persists as A^, as well as its local stable and 
unstable manifolds. 

For any general perturbation h of the form ([6]), if /'(O) = 0, the NHIM is preserved 
without any deformation for any e As = A: p = q = 0^p = q = 0. Moreover, 
if /(O) = 0, the perturbation vanishes on A, so the one-parameter family of two- 
dimensional invariant tori existing in the unperturbed case remains fixed under the 
perturbation, as in the Arnold's example of diffusion in |Arn64j . However, a generic 
perturbation, /(O) 7^ creates gaps of size ^/£ in the foliation of persisting invariant 
tori and gives rise to the large gaps problem. See section [331 for a detailed description 
of the invariant objects in this foliation. 

Although it is not a generic assumption, for the clarity of exposition and for the 
convenience of the reader not familiar with the theory of NHIM, we have chosen 
a function / in ([8]) so that the NHIM is preserved without deformation, that is 
/'(O) = 0. Nevertheless, we want to emphasize that this is not a necessary hypothesis 
for the existence of diffusion. Indeed this was not assumed in the proof of the result 
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in |DLS06l IDH09J where NHIM theory was used. Of course, any other function / 
satisfying the conditions /'(O) = and /(O) ^ wih be enough for exhibiting the 
large gaps problem, but we have chosen the concrete one ^ so that the Fourier 
coefficients for the Melnikov potential (jlSp can be computed explicitly. 

Even if for the function / in ([8]) the NHIM remains fixed, when the local stable 
and unstable manifolds are extended to global ones it is expected that in general they 
will no longer coincide and indeed they will intersect transversally along a homoclinic 
manifold. The main tool to study the splitting of the separatrix (|13|) as well as the 
associated scattering map is the Melnikov potential associated to a perturbation h 
and to the homoclinic orbit {po,qo)'- 

/H-oo 
{h{po{a), qo{cr), I,ip + Ia,s + a; 0) ,^^. 

-/i(0, 0,I,ip + la, s + a; 0))da, 
which taking into account the expression ([6|) for h, takes the form 

/oo 
[f{qo{a))-fmg{ip + Ia,s + a)da. 

In our concrete case f{q) = cosq of ([8]), the Melnikov potential turns out to be 

1 f°° 
C{I, (f,s) = - / pl{a)g{(p + Ia,s + a)da, 

^ J -oo 

and the integral can be explicitly computed by the residue theorem: 

C{I,ip,s)= ^ Ak,i{I)cos{kLp-ls- ak,i), (16) 

(fc,/)GN2 

with 

Aki{I) = 2^ ^^LzlL^-aki, (17) 

'^ ' sinhf(A;/-/) ''''' ^ ' 

where aj,^i are the general coefficients of the function g given in ^. Notice that the 
Melnikov potential ()16p has exactly the same harmonics as the perturbation (^ in ([9]). 
We now recall the role played by the Melnikov potential in the splitting of the 
separatrix ^. By Proposition 9.2 in |DLS06| . for any (/,(/?, s) G [-1/2,/*] x T^ 
and for any non-degenerate critical point r* = t*{I, ip, s) of 

r GMh^ £(/,(/? -It, s-r) (18) 

there exists a locally unique point z, 

z = z{I,ip,s;e) = {po{T*),qo{T*)J,^,s) + 0{e), (19) 

such that z£W'{As)(h VF"(Ae). 

Next, we are going to find open sets of (I, (/?, s) £ [—1/2,/^] x T^, such that the 
function ([TSD has non-degenerate critical points at r = r*(/, 93, s). 
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Taking into account the explicit expression for the Melnikov potential (J16p . the 
function (jlSp takes the form 

£(/,¥?- It, s - r) = ^ Ak,i{I)cos{k(p - Is - T{kl - I)), (20) 

(fc,0eN2 

with ^fc^i(/) as in (|17p . Notice that the Fourier coefficients Ak^i{I) are nothing else 
but the Fourier coefficients ak,i multiplied by a non-zero factor depending on kl — I 
(which decreases exponentially in \kl — l\ as \kl — l\ goes to infinity). 

The main reason for the introduction of the upper bounds for \ak,l\ in (|10P is to 
make all the computations for the series defining C{I, f, s) in (J16p and (|20p in terms 
of C^-^'{I,ip,s), its first order trigonometric polynomial in the angles {^,s). Thus, 
we have 

C{I,(p,s) = Ao^o + Ai^o{I) cos If + Aq^i cos s + 02{p,r) 

:= d^'\l,ip,s)+C>'\l,ip,s), (21) 

where Aq^ = 4ao,o! 

27r , ^ , , 27r/ 

^0.1 = ■ r^r /n\ ^o,i, and Ai,o(/) = . , . /on "i-0- ^^^ 

smh(7r/2) smh(7r/2i) 

In the formula above, without lose of genericity and to avoid cumbersome notation 
and shifts in the pictures, we have assumed that cti^q = o"o,i = 0. Otherwise, we can 
just make a translation in the variables {<p, s). 

Next we will make our computations for the function C^-^' and a posteriori we 
will justify that they are also valid for the complete function C 

So fixing {I,ip,s) S [— 1/2,I|.] x T^ we only need to study the evolution of O-^' 
along the straight lines 

R:t£R^{lP-It,s-t)£T^ (23) 

on the torus. 

By hypothesis <^TO\i . aio / and ao,i / 0, and therefore |^o,i| / and |^i^o(-^)| 7^ 
for any /. Consequently, for every fixed /, the first order trigonometric polyno- 
mial {ip,s) I—)- C^-^'{I,ip,s) possesses four non-degenerate critical points at (0,0), 
(0,7r), (vr, 0) and (vr, vr); a maximum, a minimum and two saddle points. Without 
loss of generality and for illustration purposes let us assume ai^O)'3^o,i > 0, so that 
^i_o(-^))^o,i > for any /. In this way the maximum of C^-'^'{I, •, •) is attained at 
(0,0), the minimum at (vr, vr) and the two saddles at (0, vr) and (vr, 0) (see Figure 
[T|). Of course, assuming that < p < p* and < r < r*, for p* and r* small 
enough, by the implicit function Theorem, the function C{I, •, •) possesses also four 
non-degenerate critical points close to ones of C^-^'{I, •, •), with the same properties. 

Around the two extremum points (the maximum and the minimum), its level 
curves are closed curves which fill out a basin ending at the level curve of one of 
the saddle points. Therefore, any straight line ()23p that enters into one of the two 
extremum basins is tangent to one of the closed level curves, giving rise to one or 
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Figure 1: Graph and level curves of the Melnikov potential C^-^\l,ip,s) with ai^ = 1/4, 



0-0,1 = 1/2 and 
tt/ sinh(7r/2) 



1. In this case, A 



0,0 



4, Aifl{l) = 7r/(2sinh(7r/2)) and A 



0,1 



more extrenia of (jlSp . Since the two extrenia of C{I, ■, ■) are non-degenerate for any 
/, the closed level curves close to them are convex. Therefore, any straight line ()23p 
passing close enough to the extrema gives rise to unique non-degenerate extremum of 
(jlSp . In particular, for irrational values of / the line ()23p is a dense straight line, so 
there exist an infinite number of non-degenerate extrema for (|18p . Nevertheless, and 
independently of the irrational character of /, thanks to the form of the perturbation, 
we are going to see that indeed all the closed level curves in any extremum basin are 
convex. 

To do so, we proceed in the following way. Given a fixed value of /, let us look 
for the geometric locus where the straight lines (j23p are tangent to the level curves 
of £(/,-,•). For this, we have to impose that V^^^jC, the gradient of £(/,-,•), is 
orthogonal to the slope (/, 1) of the straight line ()23p : 






(24) 



Intuitively, for fixed /, if we want to pass through a mountain of height £(/, 92, s) 
along straight lines following a direction (1, 1), equation ()24p gives the position of the 
points (99, s) of maximum height, the crest, that we denote by C = C{I). 

Using the expression ()2ip for C, equation (|24p has the form 



0. 



IAifi{I)smip + Ao,isins + 02{p,r) 
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Disregarding first the 02{p, r) term we are faced with the following implicit equa- 
tion 

q(/) sin^? + sins = 0, (25) 

where 

Mi,o(/) sinh(^/2)/^ ai,o 

Ao,i smh(7r/2/) ao,i 



Assuming that 
which holds for all / as long as 



|a(I)| < 1, 



1.03 



ain 



< 1, (27) 

ao,i 

and that p and r are small enough, equation ([2l|) defines exactly two closed curves 
Cm and Cm, parameterized by ip, which will be called crests. The crest Cm = Cm{I), 
passing through the maximum (0, 0) of C contains the saddle (0, n) and is given by 
the expression 

s = (nif, I) = - arcsin(a(/) sin ip) + 02{p, r). (28) 

The crest Cm = Cm{I) passes through the minimum (tTj-k) and contains the saddle 
(tTjO). It is given by the expression 

s = Cm{f,I) = - arcsin(a(/) sin((^ + vr)) + vr + 02{p,r). (29) 

In Figure [2] there appear these two curves (dashed black) as well as the level sets of 
the function C^-^'. 

Remark 3.1 The case when |a(/)| > 1 is analogous, but then the crests Cm and Cm 
are parameterized by the variable s. The case \a{I)\ = 1 is special, because in this 
case the union of the two curves degenerates in two straight lines along which the 
function C is constant. 

These two additional cases do not present additional difficulties for the construc- 
tion of the function r* and therefore for the subsequent mechanism of diffusion. How- 
ever, for easiness of the reading we are going to concentrate simply in the case ()27p . 
which is fulfilled thanks to the first inequality of hypothesis (|11|) . which is the only 
one that we are going to assume from now on. 

For any point (/, 99, s) G [—1/2, /^] x T^ x M, the value r* for which the function 
([T8|) has a non-degenerate critical point satisfies {I,ip — It*,s — t*) G Cm^CM- Thus, 
the non-degenerate critical points r* of (jlSp are achieved at the intersection of the 
straight line ()23p with either the crest Cm or Cm- 

Of course, for any point (/, tp, s), there may exist several intersections of the line 
([23|) with the crests Cm and Cm, giving rise to different homoclinic intersections and 
(as we will see in the next section) different scattering maps. See Figure [31 

From now on, we choose only one of these intersections, the one with the crest Cm 
and we define the function r* as the one that given (/, ip, s), associates the minimum 
\t\ such that 

{I,^-lT,S-T)eCM{I). 
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Figure 2: Closed curves satisfying (l25l) for / = 1 (the crests), dashed black, and level sets 
of the function £[-^1(1, -, ■) with the same values as in Figure [H 





Figure 3: Straight lines fl23|) with slope 1 (Left) and 0.8 (Right) and the curves fl25|) : the 
curve of the maxima (solid curve) and the curve of the minima (dashed curve). See the 
text. 
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In symbols, 

r*(/, v^, s) := min{(I, ^-lT,s-T)e Cm{I)}. (30) 

|T|eM 

To determine a domain of definition as large as possible where this function is 
continuous and to avoid a casuistic discussion, we need to check that for any fixed I, 
the straight lines cross only once the crest Cm inside the domain (ip, s) and they do 
it transversally. This implies that for any fixed /, the slope 1// of the straight lines 
is bigger than the slope of the derivative of the function ^m{^, I) for all if £ T, that 
is 

j>^{^,I), forallv^eT, (31) 

which by equations ()28p and (|26|) and assuming that p and r are small enough, is 
equivalent to 

a{I)I < 1. 

By expression (f26|) for a{I), it is easy to see that 

a{I)I < 1.6^i^, 
ao,i 

so that condition ([5T]) holds for all / as long as 

01,0 



1.6 ^^ < 1, (32) 

ao,i 

which is exactly hypotheses (jlip of Theorem 12.11 and implies ()27p . 

Under condition ()32p , one possible domain of definition H where r* is continuous 
consists of excluding, for any / E [—1/2,/^] the crest Cm{I) from the domain of 
((/9, s), that is 

H = {(/, (^, s) G [-1/2, /;] X T X M : ^(7, v?) - 27r < s < U(/, f)]- (33) 

See Figure HI 

This particular choice of r* gives rise to a homoclinic manifold to which we will 
associate in the next section an outer dynamics to the NHIM, that we will describe 
by means of the scattering map. 



Remark 3.2 Condition ([32]) is very convenient since it provides a large domain of 
definition H for t* , and therefore it allows us to define a global homoclinic manifold 
Fg. Although the condition is very restrictive, it is not necessary for the mechanism of 
diffusion. Indeed, we introduced this assumption just to avoid a casuistic description. 
However, if the condition is not satisfied we can obtain several homoclinic manifolds 
giving rise to different scattering maps, offering more possibilities for diffusion. It 
also opens the field for homoclinic bifurcations for a NHIM (see \DMR08^ ). 

Remark 3.3 The crest Cm is in the maximum basin of C, so the function C decreases 
when one travels from (0, 0) to (0, tt) increasing ip along the curve Cm o,nd increases 
when one travels from (0, vr) to (0,27r). Analogously, the other crest Cm is in the 
minimum basin of C with a complete analogous property. Since by ()24p and ()3ip . the 
curve {'^,^m{I^'^)) is never tangent to the level sets of the function C, the increase 
and decrease are strict. 
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3.2 Part 2. Outer dynamics (Scattering map) 

We have seen in the previous section that if condition ()lip is satisfied for any (I, 93, s) 
in the domain H given in (|33|) . the function (|18|) has a non-degenerate critical point 
r* given by r* = T*{I,if,s), where r* is a smooth function defined in ()30p . If < 
e < e*(I^), these critical points give rise to a homoclinic manifold F^ C VF| (ti VF^, 
along which these invariant manifolds intersect transversally. By equation ()19p . it 
has the form 

r, = {z = z{I,<f,s;e) = {po{T*),qo{T*),I,^,s)+0{e) : 
{I,ip,s)eH,T* = T*iI,^,s)£R}. 

Remark 3.4 For the experts in the splitting of separatrices, we notice that the size 
of e* required for the justification of the transversal intersection ofW? and W-' along 
Fg has to be such that the Melnikov potential ()16p gives the dominant part of the 
formula for the splitting. In our case, since C as well as its two first derivatives are 
O (exp (— 7r/2/J5l)) on the domain H , we need to impose that e* = O (exp (— 7r/2/^)) . 

In this section we will see that associated to the homoclinic manifold F^ we can 
define an outer dynamics Ss to the NHIM A and we will obtain an approximate 
explicit expression for it. 

The scattering map associated to Fg is defined in the following way: 



Ss-. H CA -^ A 

X_ I— )■ X4. 



(34) 



such that x+ = Ss{x-) if and only if there exists z G F^ such that 
dist{^t,e{z),^t,e{x±)) ^0 fort ^±00, 

where ^t,e is the fiow of Hamiltonian ([7|) . 

In words, the scattering map maps a point x_ on the NHIM to a point x+ on the 
NHIM if there exists a homoclinic orbit to the NHIM that approaches the orbit of 
X- in the past and the orbit of a;+ in the future. 

The scattering map S^ is exact and symplectic and indeed it is Hamiltonian, since 
it is given by the time e map of a Hamiltonian S^ |DLS08j . In the variables (/, (p, s) 
this implies that the following formula holds for the scattering map 



S,{I,^,s) = {I,^,s)+el—^{I,^,s),^{I,V,s),0\+O{e'). 

As it is described in equation (21) in |DH09j . the dominant term Sq of the Hamil- 
tonian Se is equal with opposite sign to the reduced Poincare function C* defined 
implicitly by 

£*(/, V3 - Is) := £(/, V3 - /r*(/, ^,s),s- t*{I, ^, s)). (35) 

It is important to notice that the Hamiltonian Sq is an autonomous 1-degree of 
freedom Hamiltonian that can be expressed in terms of the variable 9 = ip — Is: 

So{I,ip,s) = -C*{I,9), 6 = if -Is, (36) 
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Figure 4: Straight lines fl23l) with slope 1 and the curves fl25l) : the curve of the maxima 
(solid curve) and the curve of the minima (dashed curve). See the text. 



so that 

5,(/, v^, s)=[l + e-^{1, 6) + 0(^2), - e^^(/, ^) + 0{e% s J , (37) 

and the iterates under the scattering map simply follow closely the level curves of 
the reduced Poincare function ([35]) . 

In order to obtain an expression for the reduced Poincare function and for its level 
curves in our particular example, we will perform a discussion based on geometric 
considerations. 

By the definition of r*(/, (^, s) given in the previous section we have that the point 



c(/,v?,s) := {I,ip- lT*{I,ip,s),s -T*{I,ip,s)) G Cm{I), 



(38) 



belongs to the crest Cm, which is the closed curve defined in (j28|) . 

Therefore, the reduced Poincare function C* evaluated on a point (/, ip, s) in the 
domain H defined in (j33|) provides the value of the function C evaluated on c(/, (p, s), 
the closest intersection of the straight line (|'23|) starting on this point (/, f, s) with 
the curve Cm- By construction, it is clear that there is a segment of points (/, (p, s) 
in the domain H with the same c(/, (p, s) on the curve Cm- See Figure [H 

Since the function C* is constant on these segments it can be written as a function 
of only two variables: the action / and the variable 9 = ip — Is, which is 27r-periodic 
in (p and constant along the straight lines ()23p of slope l/I contained in ()33p . 

In order to obtain explicitly the expression of C* in the variable 9, we will proceed 
in the following way. Fixed /, we consider a point (/, ip* , 0) lying on the axis s = 0. 
The points (/, (p, s) in the domain H defined in (|33p satisfying 9 = (p — Is = ip* 
intersect the curve Cm on the same point c(/, ip* , 0) as (/, ip* , 0) which by ([38|) is 

{I,ip*-lT*{I,ip*,0),-T*{I,p*,0)). 
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Since each segment {9 = cte} in H intersects the s = axis at the point (1,9,0), 
the function C* has the following expression: 

C*{I,9) = C{I,9- lT*iI,9,0),-T*{I,9,0)). 

The behavior of the function C* with respect to the variable 9 (which parame- 
terizes the curve Cm) is exactly the behavior of the function C along the curve Cm, 
that we discussed in Remark 13.31 Namely, when 9 increases from to vr, one travels 
from (0,0) to (0,7r), increasing if along the curve Cm and therefore the function £.* 
decreases strictly. Equivalently, when 9 increases from vr to 27r, one travels from (0, vr) 
to (0,27r), increasing ip along the curve Cm and therefore the function C* increases 
strictly. 

We summarize some of the properties of the reduced Poincare function in the 
following Proposition: 

Proposition 3.5 For any I € [— 1/2,/ijl], the function 9 i— t- C*{I,9) has a non- 
degenerate maximum (minimum) close to 9 = (mod27r) (9 = 7r(mod27r), respec- 
tively) and is strictly monotone in all the other points 9. Moreover, it has the follow- 
ing expression 

C*{I,9) = Ao^o + Ai,o{I)cos{9 - lT*{I,9,0)) + Ao,icos{^m{I,9 - lT*{I,9,0))) 
+02{p,r), (39) 

where Ai^ and Aq^i{I) are given in ([22|) . ^m iiT' (I2SD CLnd t*{I,9,0) is defined in 
(EUD. 



Remark 3.6 Notice that the behavior of the function C* with respect to the variable 
9 is " CO sinus-like". This observation is clear when one considers the case / = 0, where 
if = 9, ^Af (0, (f) = and 

£*(0, 9) = ^0,0 + ^i,o(0) cos(^) + Ao,i + 02{p, r). 

Proposition 13.51 provides us with an exhaustive description of the level sets of the 
reduced Poincare function, giving an approximate expression in first order for the 
orbits of the scattering map S^ in ()37p . In Figure [5] we illustrate these level curves 
for a particular case. 

3.3 Part 3. Inner dynamics 

The inner dynamics is the dynamics of the flow of Hamiltonian ([7]) restricted to the 
NHIM (|12p . Indeed, by the form of the perturbation, the Hamiltonian restricted to 
NHIM ([12]) takes the explicit form 

K{I,ip,s;e) = — + eg{ip,s), (40) 

where {I,ip,s) G [-1/2, /|] x T^ and g is given in Q. 
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Figure 5: Graph and level curves of the reduced Poincare function C*{I,i 
As in Figure [T]we have chosen ai^o = 1/4 and aQi = 1/2 in (!22|) . 



given m 



In this section, we want to study the dynamics in the NHIM, that is, what are 
the invariant objects, what is the distance among them in terms of the action / and 
what are their approximate analytical expression. This section relies on the proof 
and the results of Theorem 3.1 in |DH09j . 

We already mentioned in Section 13.11 that the dynamics for the unperturbed case 
(that is, e = 0) is very simple. Indeed, all the trajectories lie on two-dimensional 
invariant tori / = const. The motion on the tori is conjugate to a rigid rotation of 
frequency vector (/, 1). Notice that the Hamiltonian K is non degenerate, that is, 






1/0. 



For £ > 0, KAM theorem ensures the preservation, with some deformation, of 
most of the invariant tori present in the unperturbed system. Indeed, those having 
frequencies "sufficiently" non-resonant, for which the influence of small divisors kl — I 
can be overcome. This is guaranteed by a Diophantine condition on the frequency 
vector: 

7 w.,. n ^ ^2\ 



\kl 



> 



miw 



V(A:,/) GZ^\{0}, 



with T > 1 and some 7 = 0{^/e) > 0. The frequency vectors (/, 1) satisfying this 
Diophantine condition fill a Cantorian set of relative measure 1 — 0{^/£), called the 
non-resonant region. 

Hence, the invariant tori with Diophantine frequencies persist under the pertur- 
bation, with some deformation. These KAM tori, which are just a continuation of 
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the ones that existed in the integrable case e = 0, are commonly known as primary 
tori and are given by the level sets of a function F of the form (see Proposition 3.24 
in [DH09]) 

F{I,^,s) = I + 0{e). 

On the contrary, the invariant tori with resonant frequencies are typically de- 
stroyed by the perturbation, creating gaps in the foliation of invariant tori of size up 
to 0{^/e) centered around resonances (indeed for the values of / such that kl — l = 0, 
for some {k, I) S N^ which is the support of the Fourier transform of the perturba- 
tion g given in Q). However, in these resonant regions, other invariant objects are 
created, like secondary tori, that is, two-dimensional invariant KAM tori contractible 
to a periodic orbit, as well as periodic orbits with associated stable and unstable 
manifolds. To prove the existence of these objects in the resonant regions and also to 
give an approximate expression for them, in |DLS061IDH09j several steps of averaging 
were performed before applying the KAM theorem to the Hamiltonian expressed in 
the averaged variables. 

More precisely, given any {ko,lo) G N^, ko / 0, gcd{kQ,lQ) = 1, for a resonant 
region centered around a resonance / = lo/ko, the invariant tori are given by the 
level sets of a function F, whose dominant term F in e is of the form (see Theorem 
3.28 in [DH09]) 

F{I,e) = M_M! +efc2f^M„Zo(^)^ (41) 

where 6 = k^ip — Iqs and the function JJ '^' " contains the resonant terms of g with 
respect to (A:o,/o)) that is, 

M 

U''°''°{e) = ^atko,tiocos{te), (42) 

t=l 

where M = e~i/(26+o) a,nd (5 is a parameter satisfying < 5 < 1/10. The lower and 
upper bounds for the coefficients aki of g provided by hypothesis ()10p ensure that the 
function (j42p reaches a non-degenerate global maximum for any / = lo/ko and any 
M > 1, independently of e. 

Indeed, the function [/'^"■'"(S) is analytic and 27r-periodic in 9. By (|10p . the 
function JJ °' " in ()42p is well approximated by its first order trigonometric polynomial 

U''°''°{e) = ak,,iocosi9) + 02{p^"r'^). (43) 

As long as p, r are small enough, and replacing 9 by 9+-K if necessary (when Uko^iQ < 0), 
the function (|43p has two non-degenerate critical points corresponding to a global 
maximum at 9 = and a global minimum a,t 9 = ir. 

Notice that in the resonant region around / = lo/ko, the angle variable 9 intro- 
duced in ()36p satisfies 9 = ip — {lo/ko)s = 9 /ko- Hence, the invariant tori in the 
resonant region can be approximated by the level sets of a function that expressed 
in the same variables (/, 9) as the scattering map in (|37p . takes the form 

F*{I,9) = ^ ^'^ + eaka,{cos{ko9) - 1), (44) 
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which is the Hamiltonian of a pendulum in the variables (/, 9 = k^O). Notice that a 
constant term — eafeg^/g has been added in order that the 0-level set corresponds to 
the separatrices of the pendulum. It is worth noticing that the size in the action I of 
the region enclosed by the two separatrices of dM]), the gap, is given by ^Je\ak~l^\. 
In terms of the variable 6, the function F* is 27rA;o-periodic and therefore the region 
enclosed by the separatrices has /cq components, the "eyes". 

The result of Theorem 3.1 in [DHOOj provides a sequence of KAM tori consisting 
of primary and secondary tori which are e^~'"''-close spaced in terms of the action 
variable I, for some r/ > 0. 

The scattering map ()34p detects the existence of heteroclinic connections among 
primary or secondary invariant tori whose distance is smaller than e. Thus, we will 
distinguish two types of resonant regions depending whether the size of the gaps 
created by the resonances is bigger or smaller than the size e of the heteroclinic 
jumps. The resonant region with big gaps corresponds to those resonances centered 
around / = Iq/Uq such that Y^i^fcoA)T ^ ^^ ^-^^ for this region we have the large 
gap problem. In these regions, we will include in the transition chain secondary tori 
given approximately by negative level sets of the function ()44p . In the regions with 
small gaps, centered around I = lo/ko such that \/\ako,io\ < ^^ , the gaps are of size 
smaller than e in terms of the I variable, so that it is possible to connect two primary 
tori on both sides of the gap. This case does not present the large gap problem and 
can be treated analogously as in the non-resonant region. Both regions, the small 
gaps region and the non-resonant region, are called in |DH09j the flat tori region 
and the dominant term in e of the invariant tori is given there in first order by the 
function 

F*iI,e) = I, (45) 

where 9 = tp — Is. 

3.4 Part 4. Combination of both dynamics 

The geometric mechanism of diffusion close to the NHIM is based on the combination 
of two types of dynamics, the inner one, provided by Hamiltonian ()40p and the outer 
one, approximately given by the time-e map of the Hamiltonian ()36p . 

Diffusion inside the three-dimensional NHIM can only take place between two- 
dimensional invariant tori. In order to overcome the obstacles of these invariant KAM 
tori present in the NHIM, we need to use the outer dynamics to "jump" from one 
KAM invariant torus to another one. Of course, for this mechanism to be successful 
the outer dynamics does not have to preserve the invariant KAM tori of the inner 
one. 

We have seen that the invariant KAM tori for Hamiltonian (j40p are given by the 
level sets of a function F, having different expressions in the flat tori region and 
in the big gaps region. Indeed, its dominant term is given by a function F* that 
depends only on the variables (1,9), having expression ()44p for the big gaps region 
and expression (|45p for the flat tori region. Moreover, in equation (|37p we showed 
that the scattering map S^ is given in first order by the time-e map of an integrable 
Hamiltonian —C*{I,9) of the form ()39p . Therefore, imposing that F* is not a first 
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integral of the Hamiltonian £*, or equivalently, that the functions F* and C* are 
functionally independent as variables of {1,0): 

dF* ■ dF*~ dF* dC* dF* dC* rr.*^*^^^ 
dl dO dl 09 de dl ^ 

guarantees that the KAM tori in the NHIM are not invariant for the scattering map. 
Moreover, it is easy to see that the Poisson bracket {F, C*} provides, in first order, 
the deformation under the scattering map of an invariant torus given by F{I, (p, s) = 
cte , since 

FoSe = F-s{F,C*} + 0{e^). (47) 

Considering the torus 

{F{I,ip,s)=E}, 

the image under the scattering map satisfies {F o S^^{I, ip, s) = E}, where F o S~^ 
has the following expression: 

FoS-^ =F + e{F,C*} + Oie^). (48) 

Using that F can be approximated by its dominant term F* , it is enough to check 

{F*,C*}^0 (49) 

to guarantee that the scattering map moves the invariant tori of the inner dynamics. 
In the flat tori region, F* is independent of e and is given in ()45p so that 



{F*,mi,e) = -^{I,d), (50) 

od 

whereas in the big gaps region F* has the form (|44p and one can easily check that 
the second term involving the derivative of F* with respect to 9 in expression ()46p 
for {F* ,C*} is small compared with the first one (see Lemma tal in |DH09j ). So, the 
dominant term in e of {F* , £*}(/, 9) is 

for any /G [-1/2,/;]. 

For Hamiltonian ([7]), by Proposition [331 we know that the reduced Poincare func- 
tion C* in (|39p is non-constant in the variable 9 and indeed, for any / G [—1/2,/;], 
the function 

~ dr.* ~ 
0^^(/,^) 

do 

vanishes only for ^ = 0, vr, and is negative for values of 9 between and vr and positive 
between n and 27r. 

Hence, in view of equation (|48p . condition ()49p guarantees that the scattering 
map (outer dynamics) moves the invariant tori of the inner dynamics. It remains to 
check now that the image under the scattering map of each of these tori intersects 
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transversally another torus of the fohation. This provides the existence of a transverse 
heterochnic connection among these tori. 

Since the different types of tori that appear in our problem have different quantita- 
tive properties and also different expressions, we consider two different cases to check 
the tranversality condition: the flat tori region and the big gaps region, introduced 
in section 13.31 

In the case of flat tori, which by equation ()45p are given approximately by / = cte , 
it is clear from expression (|48|) and (|5Up that the flat invariant tori are mapped under 
the scattering map to 

I-e^(I,0) = cte. 
89 

Because of the cosinus-like behavior of C* described in Proposition 13.51 (see also 
Remark 13. 6p , the image of this torus intersects transversally other invariant tori with 
/ = cte , that we know from the previous section that they are at a distance smaller 
than e. 

Notice that depending on the point on the torus that we choose to apply the 
scattering map, we can "jump" to different KAM invariant tori with either a higher 
or a lower value of the action /. Hence, it is possible to construct several diffusing 
orbits with increasing or decreasing /. 

In the big gaps region the computation is more complex because in this case we 
have primary and secondary KAM tori. Moreover, the invariant KAM tori are bent 
near the separatrix. In a connected component of the big gaps regions, around a 
resonance / = lo/ko, invariant KAM tori are given approximately by the implicit 
equation ^ 

F*{I,6) = E, (52) 

for F* as in (|44p and E taking discrete values around 0. 

When the value of E is positive, equation F*{I,9) = E provides two primary 
invariant tori T^ , whereas for £■ < 0, it provides a secondary torus Te- Equation 
f)52p defines two smooth surfaces given as graphs of the action / over the angle 9 
defined in a certain range: 

/ = r±{0. e) = ^ ± J2{E-eak,,i,{cos{ko9)-l)), (53) 

corresponding to the + and — sign. When £^ > these smooth surfaces correspond 
to the two primary KAM tori T^ and the graph ()53p is defined in the whole domain 
[0, 27r). For i? < 0, they correspond to the two components of the secondary KAM 
torus Te and the graph ()53p is defined in a domain strictly contained in (0, 2tt). 

By expressions (|15|) . ([^T]) and ([^5j) we have that the image under the scattering 
map of a torus satisfying {-F(/, (p, s) = E} is given in first order by the set of points 
(/, 9) satisfying 

-dC* 



F*{I,9)Te\/2iE - eak,,i,icosiko9) - 1)—^ [1,9) = E. 

d9 
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Figure 6: Invariant tori (primary and secondary) in the resonant region around / = 
(red curves) given implicitly by the level sets of the function F*{I, 6) defined in ( 144|) with 
kf) = 1, Iq = and Oi^o = 1/2- Images of these invariant tori (green curves) under the 
scattering map f l37|) generated by the reduced Poincare function C*{I,9) given in fl39|l . 



From this expression, it is clear that if the second term on the left-hand side, 

namely 

/ Fir* 

(54) 



eMl 



2{E - eako,iaicos{koe) 



— dC* 



de 



is non-constant, the image of the invariant torus intersects transversally other invari- 
ant tori at a distance smaller than e. 

In the case of primary tori (when they are defined as a graph of I over 9 in the 
whole domain [0, 27r)), there is an easy way to check this condition. Recall first that, 
by Proposition [321 the function i— t- C*{I,0) is "cosinus-like" and its dominant term 
possesses two non-degenerate critical points at = 0, vr. Therefore, one can see that 
the points on the torus corresponding to = 0, vr remain invariant for the scattering 
map, so they are not good for diffusion. However, all the other points are moved by 
the scattering map. Indeed, for 6 G (0, vr), the scattering map decreases the value of 
the action /, whereas for 6 G (7r,27r) the scattering map increases it. See Figure [6l 
Again, as in the case of the flat tori region, depending on the point on the torus that 
the scattering map is applied one can diffuse either increasing or decreasing the value 
of the action /. In this paper we are concerned with diffusion with increasing /. 

Thus, since 
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Figure 7: Illustration of how to combine the two dynamics to cross the big gaps region. 
Invariant tori for the inner dynamics (red curves) and invariant sets for the outer dynamics 
(blue curves). Inner dynamics is represented by dashed lines whereas outer dynamics is 
represented by solid lines. 



in ^ = 0, TT, because these points correspond to non-degenerate extrenia of i— )• 
£*(/, 0), it is immediate that 

^A^(^,e)/0, 
ad 

in = 0, vr. Therefore, there exists an open neighborhood around these points where 
the intersection of the invariant tori with its image under the scattering map is 
transversah This is also true for secondary tori in ^ = vr when /cq is odd. In the other 
cases and also in general, it is easy to see in an analogous way that the intersections are 
transversal because the "cosinus-like" behavior of the function Q i— t- £*(/, ^) described 
in Proposition 13.51 guarantees that 



d_ 

de 



M{i,e)^o. 



Up to this point we have shown that it is possible to construct a sequence of 
invariant tori in the inner dynamics having transverse heteroclinic connections among 
them, by means of the combination of two dynamics. See Figure [7] for an illustration 
of the combination of the two dynamics. Thus, we have constructed a transition 
chain. Finally, using a standard obstruction property (see |DLS061 IFMOlj ) one can 
show that there exists an orbit that shadows this transition chain, and Theorem 12.11 
follows. 
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